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The dynamic separation into phases of high and low baryon density in a heavy ion collision can 
enhance fluctuations of the net rapidity density of baryons compared to model expectations. We 
show how these fluctuations arise and how they can survive through freezeout. 
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QCD can exhibit a first order phase transition at high 
temperature and baryon density, culminating in a tricrit- 
ical point Specifically, below the tricritical point, 
a phase coexistence region separates distinct phases of 
QCD matter at different baryon densities, as shown in 
fig. 1. Stephanov, Rajagopal and Shuryak have pointed 
out that critical fluctuations of Et and similar meson 
measurements can lead to striking signals at the tricriti- 
cal point in relativistic heavy ion collisions 1^. 

We suggest that measurements of fluctuations of the 
net baryon number in nuclear collisions can help estab- 
lish the first order coexistence region and, ultimately, the 
tricritical point. We characterize baryon fluctuations by 
the variance Vb = (N^) - {Nb)^, where Nb = N -N is 
the net baryon number in one unit of rapidity, obtained 
from the baryon N and antibaryon N distributions; the 
average is over events. Ordinarily, net baryon fluctua- 
tions in thermal and participant nucleon models of heavy 
ion collisions satisfy 



» TV^pB/^^lB ^N + N, 



(1) 



where the second equality holds for an ideal gas ^,|^ . In 
contrast, we argue below that enhanced fluctuations oc- 
cur if the expansion of the system quenches the matter 
from an initial high density state into the phase coexis- 
tence region. At the tricritical point, these fluctuations 
diverge because VdpB/d^iB = —d'^n/dfi% ^ oo in (0), 
where is the free energy. 

Baryon fluctuation measurements add new leverage to 
a search for the flrst order region, complimenting infor- 
mation from pion and kaon interferometry, intermittency 
and wavelet analyses |^ . The latter measurements probe 
the spatial structure introduced by phase separation and 
droplet formation. By comparison, baryon fluctuations 
are weakly dependent on the morphology of the mixed 
phase, because they do not rely on distinct droplets es- 
caping a rather dense system. 

In this paper we explore the onset of baryon fluctua- 
tions and the possibility of their dissipation in nuclear 
collisions. That an order parameter such as the baryon 
density should undergo extraordinary fluctuations during 
a phase transition comes as no surprise - critical opales- 
cence results from an analogous divergence of density 
fluctuations. Furthermore, less extreme but nevertheless 
measurable density fluctuations are familiar in condensed 
matter systems that are rapidly quenched into a phase 



coexistence region 0. Perhaps more surprising is the 
possibility that observable baryon fluctuations can sur- 
vive the subsequent evolution of the system. 

As motivation, we start by describing how phase sep- 
aration can produce observable fluctuations. We then 
use a spinodal decomposition model to illustrate how 
a highly supercooled system can produce large fluctu- 
ations. Finally, we ask whether diffusion can dissipate 
these fluctuations before they are detected. Generally, 
there are two ways the dynamics can obscure large net 
baryon fluctuations in a subregion of the system. First, 
particles can diffuse throughout the fluid, diluting any 
"hot spots" and their consequent fluctuations. Second, 
fluid flow can carry the particles away to similar effect. 
We focus on the effect of diffusion because transport the- 
ory estimates indicate that the relevant diffusion coef- 
flcient can be large [^1). In principle, chemical reac- 
tions introduce dissipation by annihilating and creating 
baryon-antibaryon pairs. However, these reactions can- 
not affect the net baryon number, which is conserved, 
although they do affect the individual baryon and an- 
tibaryon fluctuations H. 




FIG. 1. Schematic QCD phase diagram ||lU2| including the 
tricritical point (p) and the spinodal region (shaded). The 
dashed curves indicate trajectories followed by ion collisions. 

Figure 1 shows the expected QCD phase diagram at 
high temperature T and baryon density pB- A phase co- 
existence region lies below the tricritical point. Inside 
this region, uniform matter must eventually break up 
into distinct domains containing the high and low density 
phases respectively. Within the coexistence region is the 
spinodal region (the shaded area). There, uniform mat- 
ter is mechanically unstable, so that small fluctuations 
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can rapidly generate bubbles throughout the fluid. Mat- 
ter is metastable between the spinodal and coexistence 
boundaries, so that fluctuations must overcome an energy 
barrier to nucleate bubbles. Phase separation by spin- 
odal decomposition and bubble nucleation have been 
discussed in the context of heavy ion collisions 
albeit with different phase diagrams in mind. 

Ideally, a heavy ion collision will produce extraordi- 
nary fluctuations if experimenters can adjust the beam 
energy and the ion combination to produce initial baryon 
densities and temperatures within the spinodal region or, 
optimally, at the tricritical point. Alternatively, fluctua- 
tions can arise if the expansion of the heavy ion system 
rapidly quenches a high density system deeply into the 
spinodal region. The dashed curves in fig. 1 show one 
such quenching trajectory (a) and one trajectory that 
only reaches the nucleation region (b). Either process 
may produce phase separation far from equilibrium . 

In spinodal decomposition, runaway density fluctu- 
ations rapidly contort and compress the high density 
fluid into filamentary regions. These runaway modes 
appear because the fluid is dynamically unstable, with 
dps/dfiB < inside the spinodal region §. These 
modes grow exponentially until the density outside the 
regions reaches the equilibrium density for that tempera- 
ture, ph, at the low-density boundary of the phase coex- 
istence curve. After this burst of nonequilbrium growth, 
phase separation in an ion collision can proceed smoothly 
as the system expands and rarefies. Nucleation may pro- 
ceed more uniformly, with perhaps one bubble growing 
smoothly as the system expands. However, there is little 
distinction between spinodal decomposition and nucle- 
ation in a highly supercooled system . 

To understand why these baryon density fluctuations 
can exceed our baseline estimate (|l|), suppose that by a 
proper time tq the phase transition has created a rel- 
atively stable initial fraction 1 — / of the low density 
bubbles within the high density phase. The net baryon 
density is then pB ~ fPq + (1 — f)ph, where pq and ph are 
the equilibrium densities of the respective phases. This 
density corresponds to a net rapidity density of baryons 
Nb ~ -^Pbtq, where A is the transverse area of the col- 
lision volume. We write 

NB~fNg + il~f)Nh, (2) 

where Nq^h = Pq,hATQ. The fluctuations of Nb are en- 
hanced relative to a uniform system because the distribu- 
tion of densities within each event is bimodal, with peaks 
at Pq and ph- The variance is therefore 

VB-Vl + f{l-f){Nq^Nh)', (3) 

where = fVq + (1 — f)Vh is the weighted average 
of the variance of each phase. If we take each compo- 
nent to be a nearly ideal gas, then Vq,h = Nqji + Nq^h, 
where N and N are the rapidity densities of baryons and 
antibaryons The quantity then reduces to (0), 



precisely the Poissonian fluctuations we would expect if 
the system were uniform. The total variance (H) exceeds 
this value by an amount proportional to the square of 
the density contrast of the two phases. A nonequilibrium 
model of the quench would be needed to compute /. Ob- 
serve that the two-phase effect, which vanishes with the 
density contrast, gives way to the critical divergence of 
(|l|) near tricritical point. 

To estimate the affect of flow on baryon evolution, we 
write the net baryon current conservation law: 

dpB/dT + d^j^ = 0. (4) 

where d/dr = u^^d^ for a fluid of four velocity u^. The 
flow of the system changes the baryon density through 
. If we neglect dissipation for the moment, then the 
net baryon number flows along with fluid as a whole, i.e. 
is ~ Pbu^- For the Bjorken scaling flow, — = 
T^^ (i, 0, 0, z) and r « \/t^ — z"^. As is well known, the 
density satisfies pb{t)t = pb{tq)tq for scaling flow. The 
rapidity density is then r independent. The variance of 
Nb for an ensemble of events i is Vb = Si(-^B ~ NbY ■ 
Differentiating, we find the scaling results 

dNB/dr w and dVs/dr w 0. (5) 

Even though the expansion after tq will cause /(t) to 
decrease, both Nb and Vb are fixed at the initial values 
It follows that (H) can represent the observed 
fluctuations, provided that the flow satisfies scaling. 

Baryon diffusion can play an important role in both 
the onset and the propagation of fluctuations. We follow 
]l3t and write jg = Pbu^ + jdiss> '^here we define so 
that the total momentum density of the fluid vanishes in 
the local rest frame. The diffusion of baryons relative to 
the fluid center of momentum gives rise to: 

jLs = -MrV OiB/T) , (6) 

where M — DdpB / dpB is loosely termed the mobility, 
D is the diffusion coefficient and = (g^'' -I- u^u'^)dfj_. 

To illustrate the onset of spinodal decomposition, we 
modify (|^) to describe the strong inhomogeneities that 
spontaneously arise in the fluid. We follow the classic lin- 
ear stability analysis of Cahn , duplicating the salient 
details here a) to motivate QCD calculations of the mi- 
croscopic inputs and h) to highlight the differences and 
similarities with DCC studies We study the evolu- 
tion of a small perturbation pB{k,t) = pB{t) ex^{ik ■ r) 
in the spinodal region, where dpB/dpB is negative. An- 
ticipating that the fluid will spontaneously become inho- 
mogeneous, we assume that the free energy of ref. [0 can 
be written as an effective Ginzburg-Landau functional of 
the baryon density, 

F{PB{r)} = J rf'V{/[pB(r1] + a^PB fm, (7) 

in the local rest frame, where / — Q/V + pbPb is the 
Helmholtz free energy density for a uniform system. Re- 
placing pLB = ipF / dNB)T.v in @ with the functional 
derivative with respect to pB, we find that 
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For dps/d^B < 0, a perturbation oi k < — 
[|5pb/9/xb|/M^]^/^ grows exponentially. The fastest 
mode at fcgp — k^j \/2 grows at the shortest time scale 



' sp 
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where /csd ~ —;=■ (9) 
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This mode dominates the early stage of the decompo- 
sition. Outside the spinodal region where OpB/dp-B > 
0, perturbations decay exponentially at rates that dif- 
fer from normal diffusion only for wavelengths smaller 
than 2n/kc (roughly the equilibrium droplet size in a 
vacuum). At these early times, the variance averaged 
over the thermal fluctuations within an event grows as 
{{pBir,t)'^)) oc e^*/'^=p. This growth exceeds that of the 
density, cx e*/'^=p, allowing large fluctuations to build. 

A quantitative treatment of spinodal decomposition 
using modern techniques Q awaits the introduction of 
kinetic terms to the coarse-grained free energy density 
iQ,^ — neither M nor ^ are known. However, we can get 
a rough idea of the time scales needed for the onset of in- 
stabilities by identifying the inverse momentum fcgp with 
the correlation length, m^^ ~ 1 fm, as in Well within 
the spinodal region, we surmise that \AI{dpB /dpB)~^\ is 
of the order of the diffusion coefficient in a stable, per- 
turbative plasma, ~ 1 — 3 fm ^. Therefore Tgp ^ 1 fm 
or smaller, suggesting that spinodal decomposition can 
erupt violently when a quench is achieved. 

One faces a similar situation in disoriented chiral con- 
densate formation There, one studies nonequihb- 
rium fluctuations of the chiral order parameter at very 
small baryon density. Equilibrium fluctuations of that 
order parameter diverge along with the baryon density 
fluctuations at the tricritical point . The essential dis- 
tinction between our first order transition and those phe- 
nomena is that, here, the low density phase can coexist 
with the high density phase in thermodynamically stable 
bubbles. This is impossible in a second order transition. 
The wavelength scale ~ 27r/fcsp of the fastest growing 
mode (^) is a nonequilibrium manifestation of the equi- 
librium size of bubbles. The second order chiral transi- 
tion is also heralded by unstable modes, but the fastest 
mode has k = 0, reflecting the critical divergence of the 
correlation length. Furthermore, the corresponding time 
scale diverges in a "critical slowing down" as the cor- 
relation length ~ oo, see, e.g., eq. (2) in [^ . 
Analogous behavior takes place at the tricritical point, 
where M~^dpB / dp,B and (|^) diverge. 

Optimistic that fluctuations can be produced by a 
quench, we turn to the dissipation of these fluctuations 
due to diffusion. For a stable liquid, we use (^,H) to ob- 
tain the diffusion equation. 



dpB/dT + pB^t,u^' = DV^PB, 



(10) 



where we take D to be constant and neglect thermodif- 
fusion. Observe that diffusion cannot affect scaling flow 



because the diffusive term vanishes if is a function of r 
alone. In that case (||) would still apply. Let us then face 
a conservative scenario in which phase separation only 
occurs as a perturbation p(r, 77) of the scaling hadronic 
system (here tanh?7 — z/t is the spatial rapidity). To be 
concrete, we take Nb and Vb to be perturbed at time tq 

by Nb ~ /(iV, - Nn) and Vb ~ /(I - f){Nq - Nh? for 
/ = /o exp{— 77^/2cr^}. Note that a « 0.88 for a spheri- 
cally symmetric perturbation. 

We describe the evolution for t > tq using ( p^ ) for 
Psd^u^ ~ Pb/t. Writing the rapidity density at fixed t 
as NB{T,rf) = pB{T,r])AT, we obtain 



T^dNB/dT = Dd^NB/dr]'^ 



(11) 



The unperturbed rapidity density Ng is constant. It fol- 
lows that Nb = NB{TQ)(l)a{il), where Nb{tq) = fo{Nq - 
Nh) and 



exp{-7?V2(o-^ + 2Ds)} 
(l-h2Dso-2)i/2 ' 



'Q 



'F ' 



(12) 



and tf is the freezeout time. Longitudinal flow lim- 
its the degree to which the Gaussian perturbation can 
be dispersed. For Tp ^ tq, we see that the rapidity 
density near ?/ « 77 « is TVs « + iVB(TQ){l + 
2i:)/(rQcr2)}-i/2 since s « Tq\ 

To study the contribution of the perturbation to the 
variance for t ^ tq, we write the event averaged Vb ~ 



V% + 2{N%Ne 



spcct to T and use (11]) to obtain 



We then differentiate Vb with re- 



T^dVB/dT 



Dd^VB/dri^ 



(13) 



in Nb/Nb- 

f2. 



We find Vb = 



to linear order 

NhfiM^i-q) - fo<l>„f^{v)}, for (f>a{v) given by (jlj). 

The diffusion coefficient for baryon current is unknown 
for the mixed phase described in |jl|,^. However, we can 
get a rough estimate for D in the pure phases from ki- 
netic theory, which implies that D ~ T^is where Vth 
is the thermal velocity of baryons and Tdifi is the relax- 
ation time for diffusion. For nucleons diffusing through a 
hadron gas, T^is ~ 35 fm and D ~ 6 fm |^ at a tempera- 
ture of 150 MeV. Flavor diffusion through a perturbative 
quark gluon plasma yields D ^ 1 — 3 fm ^ . We estimate 
diffusion in the mixed phase using the larger hadronic 
value, _D ~ 6 fm. In fact, arguments in ref. sug- 
gest similarities between hadrons and droplets of mixed 
phase. At 77 « ?/ « 0, we find the rather small decrease 
Vb/Vb(tq) - 76% for tq « 5 fm, cr = 0.88, Tp ^ 00 
and /o = 1/2. 

Experimenters can use baryon fluctuations to search 
for the tricritical point as follows. Collisions for a range 
of beam energies, ion combinations and centralities can 
produce high density systems that follow trajectories as 
in fig. 1. Results can be compared as shown in fig. 2 by 
plotting the normalized ratio llj], 
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lob^VbUN + N) 



(14) 



as a function of a normalized centrality selector, e.g. the 
total charged particle multiplicity iVch. In the absence 
of unusual fluctuations, this ratio is energy independent, 
with a value close to unity, cf. eq. (|l]). In fig. 2 we show 
the results of simulated collisions incorporating the above 
results to compute the mean and variance in the event 
generator of refs. |^,^; see these refs. for details. The 
rapidity densities of baryons, antibaryons and charged 
hadrons are taken to be 60, 15 and 300 for impact pa- 
rameter 6 = and scale with the number of participants 
for 6 > 0, as appropriate at SPS energy. A value of 
/ is assigned to each event using the ad hoc distribu- 
tion, f{b) = 0.25 [1 - [h/hof] for 6 < 6o = 3 fm. The 
mean baryon number and its fluctuations are computed 
a,t y ~ rj — Q for tq — 5 fm, Tp — 10 fm, D — 6 fm, and 
rapidity density contrasts SN — Nq — Nh — 20 and 40, 
corresponding to Pq — ph 0.10 and 0.2 fm~^ on the scale 
of normal nuclear matter density. The 'hadron' curve is 
computed assuming no enhancement. The difference be- 
tween this curve and unity is due to impact parameter 
fluctuations (volume and thermal fluctuations are 
omitted). The top curve is computed with 5N = 40 but 
without diffusion. We see that diffusion is a small effect 
for our parameter choices. 
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FIG. 2. Variance (0) as a function of charged particle 
multiplicity normalized to the average value at b = 0. 

Much more work is needed to obtain crucial infor- 
mation, such as the beam energies and rapidity ranges 
needed to probe the high density phase. To start, we 
need a better idea of the phase diagram away from the 
tricritical point. We must also understand transport pro- 
cesses in the unstable and mixed phase regions. Then, 
we can combine these ingredients to perform dynamical 
simulations such as those in Q to assess the likelihood of 
a quench. We have also overlooked the possibility of su- 
peconductivity in the high density phase, a reasonable as- 
sumption near the tricritical point Q . Superconductivity 
would modify baryon transport because a) the conden- 
sate carries baryon current and b) additional transport 
modes are possible. 

In summary, we have explored the possibility that the 
phase transition of refs. ||l|,^ can enhance fluctuations of 



the net baryon number. We find that fluctuations can 
plausibly persist through freezeout to surmount impact- 
parameter fluctuations ||l^ . A collision that reaches the 
tricritical point will have the largest fluctuations, but col- 
lisions that pass below can also be extraordinary, pro- 
vided that the dynamics can quench the system. Calcu- 
lations for y « at RHIC energy |^ suggest that proton 
fluctuations alone can reveal net baryon fluctuations at 
the level of flg. 2. Work is in progress to evaluate this 
signal at higher rapidity and lower beam energy where 
higher baryon densities are met. 
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